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This work deals with cosmological models driven by real scalar field, described by standard dy¬ 
namics in generic spherical, flat, and hyperbolic geometries. We introduce a first-order formalism, 
which shows how to relate the potential that specifies the scalar field model to Hubble’s parameter 
in a simple and direct manner. Extensions to tachyonic dynamics, and to two or more real scalar 
fields are also presented. 
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I. INTRODUCTION 

The recent discovery of cosmic acceleration Qll has deeply affected modern cosmology. One of the main 
mysteries of the cosmic acceleration is the so-called dark energy, which makes up a very significant portion of the 
total energy of the universe. The simplest way to deal with dark energy includes a cosmological constant, but we can 
also consider standard Friedmann-Robertson-Walker (FRW) models described by real scalar fields, as in quintessence 
H,.HE, 7|. For the interested reader we recommend Refs. §>, f° r some recent reviews on the subject. 

In this Letter we focus on dark energy, that is, we turn attention to FRW models described by real scalar fields. 
Our investigations inspect Einstein’s equation and the equation of motion for the scalar field in a very direct way. We 
consider models described by real scalar field in generic space-time, displaying the usual spheric, flat or hyperbolic 
spatial profile. However, we follow a very specific route, in which we use the potential of the scalar field to infer how 
the scale factor or, alternatively, Hubble’s parameter evolves in time. 

The power of the method that we develop in this Letter is related to an important simplification, which leads to 
models governed by scalar field potential of very specific form, depending on a new function, W = W {(/>), usually 
named superpotential when supersymmetry is present, which very much remind us of supergravity, although we do 
not deal with supersymmetry in the present work. As we show below, we directly relate the function W(cf> ) with 
Hubble’s parameter, and this leads to scenarios of current interest to cosmology, unveiling a new route to investigate 
the subject. Almost all calculations are done as directly as possible, and we illustrate the main results with examples 
of current interest to modern cosmology. 

Models for FRW cosmology with a single real scalar field are described by the standard action 

s = J d 4 xV=g (-\r + £(</>, d^)j (l) 

where </> describes a real scalar field and we are using 47rG = 1. The line element is ds 2 = dt 2 — a 2 (t,)dr 2 , and a{t ) is 
the scale factor. In general, the energy-momentum tensor is given by T 1 / = (p, —p, —p, —p), where p and p represent 
energy density and pressure. We use Einstein’s equation to get 

H 2 = 2p/3 — k/a 2 
d/a = — (p + 3p)/3 

where k is constant: k = 1,0, or —1, for spherical, flat, or hyperbolic geometry, respectively. 

The equation of motion for the scalar field depends on C{(j), which has the standard form 

C = ^d^d^-V{ct>) (3) 

The energy density and pressure are given by p = </> 2 /2 + V and p = </> 2 /2 — V, and the equation of motion for the 
scalar field has the form 


(2a) 

(2b) 


4> + 3Hcj) + = 0 


( 4 ) 


where H = a/a stands for Hubble’s parameter, and represents dV/dcj). 
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II. THE PROCEDURE 


We now introduce our methodology. We use the energy density, pressure and Eqs. © to get, firstly using k = 0 
for simplicity, 

H = -ft (5a) 

H' 2 = ft/ 3 + 2V/3 (5b) 

Thus, the set of Eqs. and © constitutes the equations we have to deal with in the case of scalar field with standard 
dynamics in flat geometry. 

In the standard view, since a = aft) and <fi = ftt), we have H = Hft) and from Einstein’s equation we need to see 
the potential as a function of time. However, from the equation of motion for the scalar field we have V = Vft); thus, 
to make these two views equivalent, we then need to view Hubble’s parameter as a function of the scalar field. This 
is the key point, and we make it very efficient with the introduction of a new function - W = Wft) - from which we 
can now understand that Hubble’s parameter depends on time as a function of W[(j)ft)\. That is, we write 

H = Wft) (6) 

This is a first-order differential equation for the scale factor. It allows obtaining another equation, involving the scalar 
field, in the form 

0 = -W+ (7) 

which is also a first-order equation. The two Eqs. © and © allow writing the potential as 

V = ^W 2 - (8) 

We now notice that for the above potential, solutions of the two first-order Eqs. © and 0 also solve the set of 
Eqs. © and ©. Also, the deceleration parameter is given by q = —aa/a 2 = —1 — H/H 2 , and here it has the form 
q= -1 + {W^/W) 2 . 

The above calculation directly leads to the potential in Eq. ©, which very much reminds us of supergravity, where 
supersymmetry imposes similar restriction |l2l Q10H IlfiT . However, our calculation has nothing to do with 
supergravity, which requires more sophisticated manipulations. By the way, we notice that the potential does not 
depend on the sign of W, and so the change W — * —W in the above equations in general leads to another possibility: 
H = —W and 4> = Hft. A careful search on the subject has led us to a work of Kallosh and Linde 113, in which they 
pointed out similar possibilities, although looking to the problem from another point of view, from a supersymmetric 
braneworld perspective. 

Much of the above methodology was inspired on former works, in which one uses flat space-time to investigate 
first-order equations for two or more real scalar fields 0, up to arbitrary dimensions [l9l . The positive answer 
motivates the inspection of the case with generic k. Here the set of equations is given by Eq. © and 

H = —ft + k/a 2 (9a) 

H 2 = ft/3 + 2V/3- k/a 2 (9b) 


This case is harder to implement, and we notice that the use of H = W does not suffice to solve the problem anymore. 
However, if we insist with H = W, the procedure now requires the presence of a new constraint. We follow this line, 
and we suppose that 

0 = (kaZ - Wft (10) 


where Z = Zft) is in principle arbitrary function and a is constant. This Eq. (11011 is similar to the former Eq. 0 , 
valid in the case k = 0, but now the potential has to be changed to the form 


V = -W 2 + (kaZ - Wft kaZ + -W* 


( 11 ) 


The constraint emerges consistently; it can be written as 


W^Z + W 4> Z 4l - 2 kaZZ^ - 2 WZ = 0 


( 12 ) 
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and this opens several possibilities, as we discuss below. 

We illustrate the above results with some examples. Firstly, we consider the case k = 0. We take W = A(j) n , A and 
n constants. This gives the potential 


V = ^ 2 "( 3-^) (13) 

and Hubble’s parameter is H = A 1 ^ [n(n — 2 This expression excludes the cases n = 1 and 2, which have 
to be investigated separately. The case n = 2 gives H = Ae~ 4At , and now the potential has the form V{<j>) = 
(A 2 /2)(3<)> 2 — 4)0 2 , which presents spontaneous symmetry breaking. The case n = 1 is also interesting; it reproduces 
some of the ‘negative’ potentials investigated in Ref. l2fj . 

Another example is given by W = Ae~ B ‘> > , A and B constants. This gives 

V= l -A 2 (3- B 2 )e~ 2B ^ (14) 


with H = 1/ B 2 t and q = — 1 + 1/B 2 . 

The case with k ^ 0 is more involved. We see from Eq. m that there are distinct cases to be considered. For Z = 1 
we get = 2W, and so W = Ae ± ' J ' 2 ^. This case leads to diverse decelerations, depending on the values of A,k 
and a. Another possibility is Z — . and this gives (1 — ka)= W, which shows that W is now bounded or not, 

depending on the sign of 1 — ka. The deceleration parameter in this last case is given by q— — 1 + (1 — ka)(W ( j > /W ) 2 , 
which can be lower or greater than — 1, for 1 — ka negative or positive, respectively. As an example, we choose 
W = Ae B<t> , A constant, B = ±1/Vl — ka, and a < 0. We get 

V = (1 - ka) A 2 e 2B4, (15) 

and H = 1/t. Here the deceleration parameter vanishes. 

We can also consider W = A sin (ip/y/ka — 1) for A constant and ka — 1 positive. In this case we have 

V = ^A 2 + A 2 (ka - 1) cos 2 (<t>/Vka - 1) (16) 

and H=A tanh( At). Evolution is accelerated, with q<—l. 

We introduce other examples, to further stress that the above methodology is very efficient. We consider k = 0, 
the case of flat space-time, and we suppose that W = A[B + arctan(sinh(C<£))], A, B and C constants. This gives 
the potential 


V = ‘^A 2 [B + arctan(sinh(C<(>))] 2 — -A 2 C 2 sech 2 (C(/)) 


(17) 


We plot this potential for some values of A, B , and C in Fig. [1], We illustrate some possibilities with A = 1, C = -y/2 
and: B = 0, which gives a symmetric potential, similar to the form required in the bicycling scenario of Ref. |20l : 
B £ (0, 7t/2) which gives asymmetric potentials, similar to the case required for cyclic evolution B = 7r/2, 

which gives a kink-like potential, which drives the model to the case of static universe. 

We now investigate Hubble’s parameter for the above model. The calculation is easy, direct; it allows writing 
H = A[B — arctan (AC 2 t)], which shows that A, B and C control the way H evolves in time. 

Another example in flat space-time is given by [lflj 


w * = c '-2f=-i+' + wri l - 


(18) 


where C p = 2p/{Ap 2 — 1) and p = 1,3,..., is odd integer. The potential is awkward, but now Hubble’s parameter is 
given by 

H p = Cp~ P tanh 2p ~ 1 {t/p) + P tanh 2p+1 (f/p) (19) 

2p — 1 2p + 1 

which is plotted in Fig. [2] for some values of p. 
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FIG. 1: The potential of Eq. 11 71 for A = 1, C = V2, and B = 0,7t/4, and 7t/2, plotted with thin, thick, and thicker lines, 
respectively. 



FIG. 2: Hubble’s parameter of Eq. m for p = 1, 3, and 5, plotted with thin, thick, and thicker lines, respectively. 

III. TACHYONIC DYNAMICS 

The power of the method can be further extended to the case of scalar field with tachyonic dynamics; see for 
instance Refs. for specific details on scalar fields with tachyonic dynamics, related to the subject under 

investigation. The Lagrange density has the form 


£ t = -!/(</>) V l-d^d^ (20) 

The energy density and pressure are given by p t = V/(l — if) 2 ) 1 / 2 and p t = —V(l — (j) 2 ) 1 / 2 . We consider k = 0, and 
we use Eqs. 0 , the energy density and pressure, and the equation of motion for the scalar held to get 

<j> + (1 — <p 2 )(3H<j) + V^/V) = 0 
H = -V<j) 2 /{1 - <fi 2 ) 1/2 
H 2 = 2E/3(1 - </> 2 ) 1/2 


(21a) 

(21b) 

(21c) 
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The case of tachyonic dynamics is more involved, but we keep using H = W. This leads to the equation 


</> = 


2 

3 W 2 


( 22 ) 


and the potential is now given by 


Here the deceleration parameter is q = — 1 + 2W 2 /3W A . 

We consider an example, in which W = A/cfi. We use Eq. (12311 to obtain 


(23) 




(24) 


which requires that A > 2/3. We use this result to obtain H = 3A 2 /2t, and now q = — 1 + 2/3 A 2 . This example 
displays the inverse square potential, which is required for the presence of scaling solutions; see [25| and references 
therein for other details. 


IV. TWO-FIELD MODELS 

We now extend the above procedure to the case of two or more real scalar Eelds with standard dynamics. We 
consider two-field models for simplicity, and we investigate the important case of flat geometry. Here we have to 
change the Lagrange density in Eq. 0 to the form 

£2 = - V (</>, x) (25) 

and now we get the new set of equations 

4> + 3 Hi + = 0 

X + 3 Hx + V x =0 
H = —4> 2 - X 2 
H 2 =<p 2 /3 + x 2 /3 + TV 13 

As before, we insist with H = W, but now W = W((f>,x) suggests that we write the two first-order equations 

</> = -W+, X = -W x (27) 

and the potential 

V = lw 2 - l -W 2 - ±W 2 (28) 

It is not hard to show that solutions of the above first-order equations also solve the set of Eqs. Also, here the 
deceleration parameter is q = — 1 + ( W^/W ) 2 + ( W x /W) 2 . 

The above procedure can be generalized to three or more real scalar fields straightforwardly, and it opens interesting 
possibilities for coupled fields in flat geometry. For instance, if we consider an additive W. that is, if we take 
^(^X) = + W 2 ( X ), we get the potential in the form V(<j>,x) = + V 2 (x) + 3Wi(^)W2(x)- It shows that 

the interactions appear as the product of the two independent Wi and W 2 . Here we see that Hubble’s parameter is 
also additive, that is, H = H\ + H 2 . However, acceleration is not additive anymore. 

We illustrate the case of two fields firstly with W which involves product of the fields. An example which we solve 
exactly is given by W = A(f> + Bx — C(j>x, for A, B , and C constants. The potential is 

V = -±(A 2 + B 2 ) + C(A X + B<?>) + ^(3A 2 -C 2 )c/> 2 

+ \(3B 2 - C 2 )x 2 - 3C{A<!> + B X )4>X + \c 2 <t > 2 X 2 


(26a) 

(26b) 

(26c) 

(26d) 
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The first-order equations are <f> = —A + C\ and \ + C*<(>. The solutions require two new constants, D and E, 

and they give H = AB/C + CE 1 2 3 4 5 e~ 2Ct - CD 2 e 2dt . 

Another example is for W additive, given by W = Acj) 2 + W p (y), A constant, and W p as in Eq. (ITsl) . The potential 
is awkward, but now Hubble’s parameter is 


H = Ae~ AM + H p (29) 

where H p is given by Eq. iPJ) . We plot this parameter in Fig. [3] to illustrate its behavior in terms of p , for A = 1. 
A comparison between Fig. [2] and Fig. [3] shows one field affecting the other behavior explicitly. 

Evidently, we can use other models with several fields to get to richer scenarios, as we will explore in another work. 
Here we learn from the above examples how the scalar fields directly contribute to modify the cosmic evolution, and 
this is known to be of interest to modern cosmology |27| ■ 



FIG. 3: Hubble’s parameter of Eq. EM for A = 1 and p = 1, 3, and 5, plotted with thin, thick, and thicker lines, respectively. 


V. ENDING COMMENTS 

In this work we have shown how to write a first-order formalism to FRW cosmology, described by a single real scalar 
field with either standard or tachyonic dynamics for generic spherical, flat or hyperbolic spatial geometry. The crucial 
ingredient was the introduction of a new function, W — W((/)), from which we could express Hubble’s parameter in 
the form H(t) = W[(f>(t)\. The importance of the procedure is related not only to the improvement of the precess of 
finding explicit solution, but also to the opening of another route, in which we can very fast and directly write H(t ) 
once W (</>) is given. As we have shown, the present investigations are of direct interest to modern cosmology, since 
they seem to open several distinct possibilities of investigation. 

The interest in the subject broadens with the extension of the method to the case of several fields, since now we can 
find models in which one field can be used to affect the behavior of other fields, unveiling the possibility to control a 
given phase and to link different phases of the cosmic evolution. 

The authors would like to thank Francisco Brito and Carlos Pires for discussions, and CAPES, CNPq, 
PADCT/CNPq, and PRONEX/CNPq/FAPESQ for partial support. 
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